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SUMMARY 


A  number  of  the  most  promising  integral  methods  for  solving  approxi¬ 
mately  the  compressible -laminar-boundary -layer  equations  are  investigated 
in  order  to  determine  a  computationally  convenient  and  sufficiently 
accurate  method  of  calculating  boundary-layer  characteristics.  The  chief 
methods  considered  are:  (a)  The  one-parameter  Karman-Pohlhausen  method, 
with  three  different  assumptions  for  the  velocity  profiles,  and  (b)  the 
two -parameter  method,  first  applied  by  Sutton,  with  two,  different  assump¬ 
tions  for  the  velocity  profiles.  After  the  methods  are  explicitly 
described  in  general  terms  for  the  case  of  zero  pressure  gradient  and 
for  the  case  of  a  pressure  gradient  in  the  direction  of  flow  with  zero 
heat  transfer,  they  are  applied  to  the  calculation  of  the  compressible 
laminar  boundary  layer  over  a  surface  with  zero  pressure  gradient,  with 
and  without  heat  transfer  at  the  surface,  for  the  purpose  of  establishing 
the  accuracy  of  the  methods.  Comparison  of  the  results  is  made  with  those 
of  known  exact  solutions  for  skin-friction  and  heat-transfer  coefficients, 
velocity  profiles,  velocity  derivatives,  and  especially  laminar-boundary- 
layer  stability.  From  this  comparison  it  is  found  that  the  Karman- 
Pohlhausen  method  with  a  sixth-degree  polynomial  as  the  velocity  profile 
is  the  most  suitable  for  many  practical  purposes. 


INTRODUCTION 


It  is  well-known  that  the  differential  equations  of  two-dimensional 
compressible -laminar -boundary-layer  flow  are  difficult  to  solve  exactly. 
Stewartson  (reference  1)  and  Illingworth  (reference  2)  have  recently  shown 
that  if  the  Prandtl  number  is  unity  and  the  viscosity  coefficient  is 
proportional  to  the  temperature,  then  the  equations  for  the  compressible 
he at -insulated  boundary  layer  with  a  given  pressure  gradient  can  be  trans¬ 
formed  into  the  equations  for  an  incompressible  boundary  layer  with  a 
different  pressure  gradient}  however,  this  principle  appears  at  present 
tedious  to  apply  in  practice. 


The  most  frequently  used  and  most  fruitful  methods  of  solving  the 
boundary-layer  equations  approximately  are  the  integral  methods,  in 
which  the  partial  differential  equations  are  integrated  over  the  boundary- 
layer  thickness,  and  are  hence  satisfied  only  "in  the  average."  By 
assuming  definite  forms  for  the  velocity  profiles  as  functions  of  the 
normal  distance,  ordinary  differential  equations  are  obtained,  with 
distance  along  the  surface  as  the  independent  variable.  Any  integral 
method  may  be  regarded  as  either  of  two  types:  (a)  The  single -integral 
type ,  in  which  the  partial  differential  equations  are  integrated  once 
across  the  boundary-layer  thickness,  and  the  profiles  contain  a  single 
parameter  to  be  determined  by  the  resulting  ordinary  differential  equa¬ 
tion}  (b)  the  multiple -integral  type,  in  which  several  (say  m)  integral 
equations  are  used,  and  the  assumed  velocity  profiles  contain  m  param¬ 
eters  to  be  determined  by  the  m  resulting  ordinary  differential  equations. 

The  best-known  integral  method  is  the  Karman-Pohlhausen,  which  is 
of  type  (a)  with  fourth-degree  profiles.  This  method  has  been  found 
quite  useful  for  incompressible  flow  (cf.,  e.g.,  Dryden,  reference  3) 
but  it  has  two  important  limitations.  It  fails  to  predict  the  separation 
point  accurately  in  an  adverse  pressure  gradient  and  it  often  does  not 
give  sufficiently  accurate  ‘results  for  stability  calculations  based  on 
criteria  developed  by  Lin  and  Lees  (references  k  and  5).  For  these 
purposes  refinements  in  the  usual  Karman-Pohlhausen  method  must  be  made. 

In  the  methods  of  type  (a)  the  refinements  usually  consist  of 
assuming  types  of  profiles  which  satisfy  more  boundary  conditions  than 
the  fourth-degree  profiles.  Schlichting  and  Ulrich  (reference  7),  for 
example,  have  used  sixth-degree  profiles  for  incompressible  flow,  satis¬ 
fying  an  additional  boundary  condition  at  the  wall  and  also  at  the  outer 
edge  of  the  boundary  layer.  Satisfactory  results  were  obtained,  except 
for  flow  in  the  vicinity  of  a  .stagnation  point.  Weil  (reference  8)  has 
recently  applied  this  method  to  compressible  flow  with  zero  heat  transfer. 
However,  no  investigation  was  made  here  of  the  expected  accuracy  of  the 
results,  although  it  was  pointed  out  by  Weil  that  the  use  of  sixth-degree 
profiles  is  expected  to  yield  satisfactory  results  for  stability  calcu¬ 
lations  since  these  calculations  involve  first  and  second  derivatives  of 
the  velocity,  and  the  velocity  profiles  are  made  to  satisfy  additional 
conditions  involving  the  rate  of  change  of  second  derivatives  (viz,  third 
derivatives) . 

Timman  (reference  9)  has  suggested  the  use  of  exponential  profiles 
based  on  exact  profiles  derived  from  the  solutions  of  Von  Karman  and 


^Lees  (reference  6)  has  recently  applied  the  Karman-Pohlhausen 
method  to  an  investigation  of  the  stability  of  compressible  laminar 
boundary  layers  with  favorable  pressure  gradients,  but  this  investigation 
was  considered  by  its  author  to  be  essentially  qualitative. 


NACA  TN  2655 


3 


Millikan  (reference  10).  Although  the  calculations  are  thereby  made 
more  tedious,  it  is  believed  that  no  significant  improvement  over  poly¬ 
nomial  profiles  should  be  expected,  since  the  latter  type  of  profile 
usually  satisfies  a  fairly  large  number  of  conditions  at  the  outer 
boundary-layer  edge  anyway,  and  since  primary  interest  usually  lies  in 
the  region  near  the  wall.  Moreover,  Yuan  (reference  11)  and  Lew  (refer¬ 
ence  12)  have  also  used  exponential  profiles  with  no  evident  improve¬ 
ment  in  accuracy.  In  cases  of  adverse  pressure  gradients,  Timman  has 
suggested  the  use  of  a  special  profile  satisfying  an  additional  condi¬ 
tion  involving  the  fourth  derivative  of  velocity  at  the  separation  point. 
Considerable  improvement  in  the  accuracy  of  prediction  of  the  separation 
point  was  thereby  obtained  (reference  9). 

Loitsianskii  (reference  13)  has  suggested  a  modification  of  the 
Karman-Pohlhausen  method  based  on  multiplying  the  momentum  equation  by 
a  small  variation  of  velocity  and  then  integrating  across  the  boundary- 
layer  thickness.  A  velocity  profile  with  a  single  undetermined  param¬ 
eter  is,  as  usual,  assumed.  The  skin  friction  is  subsequently  calculated 
by  means  of  the  Karman  momentum  equation.  The  method  was  applied  in 
reference  13  for  several  cases  of  incompressible  flow,  but  the  results 
did  not  seem  to  indicate  superiority  of  this  method  over  the  usual  Karman- 
Pohlhausen  method  with  a  fourth-degree  profile. 

In  the  integral  methods  of  type  (b)  the  Karman-Pohlhausen  method 
is  extended  by  deriving  more  than  one  integral  equation.  This  can  be 
done  by  multiplying  the  momentum  partial  differential  equation  by  a 
series  of  different  factors,  and  then  by  integrating  the  resulting  equa¬ 
tions  over  the  boundary -layer  thickness.  The  factors  which  have  usually 
been  chosen  are  integral  powers  of  either  the  velocity  (Leibenson,  refer¬ 
ence  Ik)  Golubev,  reference  15;  Sutton,  reference  l6;  and  Wieghardt, 
reference  17)  or  the  normal  distance  (Whitehead,  reference  18).  A  further 
possibility,  suggested  and  applied  by  Whitehead,  is  successive  integra¬ 
tion  of  the  momentum  equation.  The  integral  equations  obtained  by  these 
procedures  are  to  some  extent  analogous  to  those  which  would  be  obtained 
by  the  method  of  moments,  and  an  infinite  number  of  such  equations  would 
be  equivalent  to  the  original  partial  differential  equation.  Because  of 
the  elaborate  nature  of  the  calculations  required  in  such  a  procedure, 
however,  only  the  first  two  of  the  infinite  set  of  equations  have  usually 
been  considered. 

The  multiple -integral  methods  have  thus  far  been  developed  and 
applied  only  for  incompressible  flow.  In  the  present  investigation  only 
the  use  of  powers  of  velocity  as  factors  will  be  considered.  In  the 
application  made  by  Sutton  (reference  l6)  a  fourth-degree  velocity  pro¬ 
file  was  assumed  with  two  undetermined  parameters.  However,  one  of  the 
boundary  conditions  at  the  wall  ordinarily  satisfied  in  the  Karman- 
Pohlhausen  method  was  not  satisfied.  Wieghardt  (reference  17)  has  also 
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used  this  two -parameter  method,  but  he  has  assumed  eleventh-degree 
velocity  profiles  satisfying  additional  conditions  at  the  wall  (in  fact, 
the  same  as  those  satisfied  by  the  sixth-degree  profiles  of  references  7 
and  8)  and  at  the  outer  boundary-layer  edge.  Although  results  of  com¬ 
paratively  high  accuracy  can  thereby  be  obtained,  Wieghardt’s  method 
can  be  quite  tedious  in  practice.  This,  in  fact,  is  one  of  the  general 
disadvantages  of  the  multiple -integral  methods. 

The  general  case  of  heat  transfer  in  a  compressible  laminar  boundary 
layer  with  a  pressure  gradient  is  complicated  by  the  fact  that  in  this 
case  there  is  no  known  solution  of  the  energy  partial  differential  equa¬ 
tion  giving  the  temperature  explicitly  as  a  function  of  the  velocity. 

One  means  of  treating  such  a  case  is  by  transforming  the  energy,  as  well 
as  the  momentum,  equation  into  a  differential-integral  equation,  and 
assuming  a  profile  not  only  for  the  velocity  but  also  for  the  temperature. 
Kalikhman  (reference  19)  has  investigated  this  case  by  this  means  using, 
analogously  to  the  ordinary  Karman-Pohlhausen  method,  fourth-degree  pro¬ 
files  for  both  velocity  and  stagnation  enthalpy.  Although  important 
useful  results  are  thus  obtained,  their  accuracy  is  subject  to  the  limi¬ 
tations  of  the  Karman-Pohlhausen  method  previously  discussed. 

The  aim  of  the  present  study  is  to  investigate  the  practical  feasi¬ 
bility  of  the  most  promising  integral  methods  (single -integral  and  double¬ 
integral  methods  with  various  polynomial  velocity  profiles)  from  the  point 
of  view  of  simultaneous  accuracy  and  ease  of  computation.  The  approach 
here  is  primarily  a  posteriori.  The  implications  of  each  of  the  methods 
considered  are  developed  for  only  the  simplest  case,  that  of  a  surface 
with  zero  axial  pressure  gradient  in  a  subsonic  and  in  a  supersonic 
stream,  and  these  are  compared  with  the  corresponding  implications  of  an 
exact  solution.  In  particular,  results  for  skin  friction  and  heat  trans¬ 
fer  at  the  surface,  velocity  profiles,  velocity  derivatives,  and  laminar 
stability  based  on  the  work  of  Lees  (reference  5)  have  been  considered. 

The  comparison  of  stability  is  particularly  critical,  since  it  is  here 
that  the  largest  errors  are  incurred  in  the  approximate  methods,  and 
that  the  greatest  differences  among  the  results  of  the  various  methods 
appear. 

The  following  points  are  among  the  distinguishing  features  of  this 
study:  (a)  The  two-parameter  method  involving  two  differential-integral 

equations  is  developed  and  applied  to  compressible  flow,  and  (b)  the 
implications  of  the  various  methods  are  directly  compared  with  respect 
to  stability  criteria;  this,  moreover,  is  the  main  basis  here  for  judging 
the  relative  accuracy  of  the  methods  for  practical  purposes. 

Although  it  does  not  necessarily  follow,  without  further  investi¬ 
gation,  that  the  conclusions  drawn  from  the  analysis  of  the  solutions 
for  the  flow  with  zero  axial  pressure  gradient  will  be  valid,  without 
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modification,  for  flows  with  axial  pressure  gradients,  it  is  believed 
that  the  comparisons  and  resulting  conclusions  should  nevertheless  serve 
as  an  indication  of  the  relative  merits  of  the  various  methods  in  more 
general  cases. 

It  may  be  noted  that  Mangier* s  tranformation  (reference  20),  which 
brings  the  boundary-layer  equations  of  axially  symmetric  flow  into  the 
form  of  the  two-dimensional  equations,  further  extends  the  usefulness 
of  two -dimensional -flow ‘solutions.  Moreover,  axially  symmetric  flows 
of  constant  pressure  thereby  lead  to  equations  analogous  to  those  of 
two-dimensional  flow  with  zero  pressure  gradient. 

This  work,  carried  out  at  the  Polytechnic  Institute  of  Brooklyn 
Aeronautical  Laboratories,  was  sponsored  by  and  conducted  with  the 
financial  assistance  of  the  National  Advisory  Committee  for  Aeronautics. 
The  authors  wish  to  acknowledge  the  helpful  discussions  of  this  research 
with  Professor  R.  Paul  Harrington  and  the  valuable  contribution  of 
Dr.  H.  G.  Lew  to  the  formulation  of  the  program  in  its  early  stages. 


SYMBOLS 


aj 


h 

k 

M 

%u 


coefficients  appearing  in  velocity  profile 
factor  of  proportionality  in  equation  iL  =  C— 

^00 

coefficient  of  average  skin-friction  drag  for  surface  of 
length  x 

specific  heat  at  constant  pressure 
specific  heat  at  constant  volume 
integrals  defined  in  equations  (22) 

local  heat-transfer  coefficient 
thermal  conductivity  of  fluid 
Mach  number 

local  Nusselt  number  fhx/k  ) 
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N 


Pr 


cr 


Re 

s 

t 

T 

U1 

Uco 

U,V 


€ 

P 

P 

7 

V 

6 

5t 


Prandtl  number  of  gas  ^jiCp/kj 
pressure 

gas  constant  -  cv^ 

critical  Reynolds  number 

Reynolds  number  based  on  length  x  J^Vw| 

Reynolds  number  based  on  momentum  thickness  0 

Sutherland1 s  constant 
transformation  variable  defined  in  equation  (15) 
absolute  temperature 

velocity  at  outer  edge  of  boundary  layer 
free -stream  velocity 


velocity  components  in  x-  and  y-directions,  respectively 

coordinates  parallel  and  perpendicular  to  surface, 
respectively 

variable  defined  in  equation  (35) 
mass  density 

absolute  coefficient  of  viscosity 
ratio  of  specific  heats  ^Cp^c^ 

kinematic  viscosity  (n/p) 

thickness  of  boundary -layer  velocity  profile  in  xy -plane 
thickness  of  boundary-layer  velocity  profile  in  xt-plane 


w 


* 


* 


T 


transformation  variable 


NACA  TN  2655 


7 


Subscripts: 

00  free-stream  conditions  for  subsonic  flowj  uniform  conditions 

behind  shock  for  supersonic  flow 

1  local  conditions  at  outer  edge  of  boundary  layer 

o  conditions  at  surface 

c  values  at  critical  coordinate  for  stability 

BASIC  EQUATIONS 

The  following  equations  describe  the  steady,  two-dimensional,  laminar- 
boundary-layer  flow  of  a  compressible  gas  along  a  surface  whose  radius 
of  curvature  is  large  compared  with  the  boundary -layer  thickness: 

Momentum  equation  in  x-direction: 

pu(du/dx)  +  pv(du/dy)  =  -(dp/dx)  +  ( (l) 
Momentum  equation  in  y-direction: 

dp/dy  =0  (2) 

Equation  of  continuity: 

(d/ck)(pu)  +  (d/dyUpv)  =  0  (3) 

Equation  of  state: 


p  =  pRT 


Energy  equation: 

pucpOT/^x)  +  pvcp(dT/dy)  =  u(dp/Bx)  +  (d/cfy)[k(dT/5y)]  +  p(du/dy)2  (5) 
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As  a  consequence  of  equation  (2)  the  static  pressure  p  at  any 
axial  position  x  in  the  boundary  layer  is  equal  to  the  corresponding 
static  pressure  in  the  potential  flow,  in  which  conditions  are  assumed 
to  be  known. 

From  the  momentum  (Bernoulli)  equation  in  potential  flow,  the 
pressure  p  can  be  expressed  in  terms  of  the  potential-flow  velocity 
distribution,  thus: 


dp/dx  =  dp/dx 

=  -PlUlV  #  (6) 


where  the  prime  denotes  differentiation  with  respect  to  x,  and  the  sub¬ 
script  1,  conditions  at  the  outer  edge  of  the  boundary  layer.  As  is 
usual  in  aeronautical  problems,  the  coefficients  of  viscosity  and  heat 
conductivity  are  taken  to  be  known  functions  of  the  temperature  only, 
while  the  coefficient  of  specific  heat  at  constant  pressure  will  be 
assumed  constant. 

The  treatment  of  these  equations  is  greatly  simplified  if  use  is 
made  of  well-known  algebraic  velocity- temperature  relations  which  are 
exactly  valid  for  a  Prandtl  number  of  unity.  These  relations  are  for 
the  case  of  zero  heat  transfer  with  or  without  an  axial  pressure  gradient 


+  cpT  =  A 


(7) 


and  for  the  case  of  zero  axial  pressure  gradient  with  heat  transfer 


+  CpT  =  A  +  Bu 


(8) 


where  A  and  B  are  arbitrary  constants  determined  from  wall  and  free¬ 
st  re  am  conditions.  With  the  use  of  these  conditions,  equation  (7)  becomes 


t/t-l  =  1  +  (7  -  1) 


(Ml2/2)[l  -  (u/Ui)2 


m 


(9) 
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while  equation.  (8)  becomes 

i/*i  -  (yTi)  -|(to/ti)  -  [L  +  (r  -  i)(mi2/||  (u/uiXr 

(7  -  1)  (m12/2)(u/u1)2  .  (10) 

The  use  of  equations  (9)  and  (10)  which  are  valid  only  for  a  Prandtl 
number  of  1  in  lieu  of  the  more  general  energy  equation  (5)  is  justified 
primarily  by  the  simplifications  introduced  in  the  analysis.  However, 
it  has  been  found  that  momentum  boundary -layer  characteristics  such  as 
skin  friction  are  not  significantly  affected  by  a  change  in  Prandtl 
number  from  unity  to  the  actual  value  for  air,  which  is  between  0.65 
and  O.76.  Heat-transfer  coefficients  are  affected  approximately  10  per¬ 
cent  by  such  a  change. 

Before  discussing  the  approximate  solution  of  the  basic  equations 
by  means  of  integral  methods,  it  is  convenient  to  present  several  sub¬ 
sidiary  relations  for  later  use.  Since  the  pressure  within  the  boundary 
layer  at  any  given  value  of  x  is  a  constant,  equation  (4)  gives 


p/pi  =  Ti/T 


(id 


The  most  accurate  representation  of  the  variation  of  the  coefficient  of 
viscosity  of  gases  with  temperature  is  usually  considered  to  be  Sutherland* s 
formula,  namely, 

^00=  (T/Tco)3/2(Tc,+  S)AT+S)  <12) 

where  S  is  a  constant  which  for  air  is  216°  R  and  where  the  sub¬ 
script  00  denotes  reference  conditions  taken  in  the  undisturbed  free 
stream  for  subsonic  flow  and  in  the  uniform  flow  behind  any  shock  wave 
associated  with  the  leading  edge  for  supersonic  flow.  Although  it  is 
possible  in  principle  to  use  equation  (12)  in  the  methods  employed  here, 
it  is  inconvenient  to  do  so  since  certain  integrals  must  then  be'  evaluated 
numerically.  It  is  thus  desirable  that  some  acceptable  simplifying 
approximation  to  equation  (12)  be  made.  In  reference  (21)  Chapman  and 
Rubesin  suggested  that  the  mathematical  advantages  of  assuming  p  and 
T  to  be  linearly  related  can  be  obtained  along  with  satisfactory  accuracy 
if  it  is  assumed  that 
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•*M.  -  C(T/T„) 


(13) 

V 


where  C  is  a  factor  which  is  chosen  so  that  the  Sutherland  viscosity- 
temperature  relation  is  exactly  satisfied  at  the  wall  temperature  T0. 
Thus, 


C  =  (To/TJl/2(Tro  +  s)/(tq  +  S) 


(1^) 


Equation  (13)  will  be  used  throughout  this  report. 

It  is  convenient  in  the  method  used  here  to  apply  the  Dorodnitzyn 
transformation.  Thus  a  new  variable  t  is  defined  so  that  for  a  given 
value  of  x  the  physical  variable  y  is  given  by  the  equation 


y 


dt 


(15) 


Correspondingly  the  thickness  5-j-  of  the  boundary  layer  in  the  xt-plane 
is  defined  as  the  value  of  t  when  y  =  5,  and  therefore 


6 


(16) 


With  the  additional  relations  given  by  equations  (11),  (13), 

(14),  and  (15),  equation  (l)  can  be  recast  into  an  infinite  set  of 
integral-differential  equations.  This  can  be  done  by  multiplying  equa¬ 
tion  (l)  by  u11  dy(n  =  1,2,3,  .  .  .)  and  by  integrating  from  y  =  0  to 
y  =  B  with  the  boundary  conditions  that  at  y  =  0 


and  at  y  =  B 


u  =  Ux 


(17) 


6u/c)y  =  0 
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After  transforming  to  the  xt -plane,  nondimensionalizing,  and  introducing 
T  =  the  following  set  of  equations  is  obtained: 


(n  +  l)p1U1 


n+2 


-1 


Bt(d/dx) 


piuin+2Bt /  Wul)  _ 


1 


-  1) 


n+1 1  dT^ 


KHl'tW /J"/^  -  (,T/T i)(u/ni)n 


dT  = 


C^p^T^”1  Ra/Br)  (u/u^ 


dT  n  >  1 


n  =  0 


(18) 


The  application  of  the  infinite  set  of  equations  represented  by 
equations  (l8)  will  now  he  discussed.  By  means  of  equations  (9)  and  (10) 
the  temperature  can  he  expressed  in  terms  Of  the  velocity  component  u. 
Now  if  it  were. assumed  that 


00 

u/Ui  = 

i=o 


then  several  of  the  a  j  coefficients 

the  other  coefficients  so  that  certain  boundary  conditions,  at  least 
those  in  the  xt-plane  corresponding  to  those  given  hy  equations  (17)  in 
the  xy -plane ,  would  he  satisfied.  The  rest  of  the  coefficients  and  6^ 

could  he  calculated  from  the  infinite  set  of  first-order  ordinary  differ¬ 
ential  equations  represented  hy  equations  (l8).  Such  a  solution  would 
in  principle  he  exact  and  would  he  equivalent  to  a  solution  of  the 
original  partial  differential  equations. 

In  actuality,  however,  since  the  attainment  of  such  a  solution 
would  he  quite  difficult,  it  will  usually  he  considered  sufficient  for 
practical  purposes  to  obtain  approximate  solutions  to  the  infinite  set 
hy  taking  only  the  first  few  differential  equations  of  this  set  and 
correspondingly  only  the  first  few  terms  of  the  power  series  in  equa¬ 
tion  (19).  Indeed,  much  information  of  engineering  importance  for  the 


aj(x)T  ^  (19) 

could  he  determined  in  terms  of 


4 
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flat-plate  case  has  been  obtained  in  the  past  by  using  only  the  equation 
corresponding  to  n  =  0,  that  is,  to  the  Von  Karman  momentum  integral 
(e.g,,  references  3>  11 >  and  12)*  Additional  accuracy  may  be  obtainable 
by  the  use  of  a  second  equation,  that  is,  by  the  use  of  equations  (l8) 
for  n  =  0  and  n  =  1.  Therefore,  these  two  equations  will  now  be 
explicitly  written. 

For  n  =  0: 

(Fl/2)(6t2)'  +  [V  +  Fl(l°Se  Pl»l2)'  -  f2(loee  Ui)]^2  = 
(F3|i.CTl)/(PlV«.)  (2°) 

and  for  n  =  1: 

(Flt/‘*)(6t2)  +  [(V/2)  +  (FW2)  (l0«c  Plul3)'  -  f5  (Foge  ul)]  \2  - 

[(W)/(Pl»L)](Tl/*.)  (a) 

where  the  prime  denotes  differentiation  with  respect  to  x  and  where 
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In  these  Fj  integrals  either  equation  (9)  or  (lO)‘,  as  the  case  may 
he,  is  used  to  express  T/T^. 

In  the  actual  application  of  either  or  both  equations  (20)  and  (21), 
u/u is  assumed  as  some  function  of  t.  For  convenience  a  polynomial 

in  t  is  here  assumed,  that  is 


where  N  is  the  degree  of  the  assumed  polynomial  containing  (N  +  1)  of 
the  aj's. 

The  ajfs  are  determined  from  the  boundary  conditions  and  from  the 

differential  equations  (20)  and  (2l).  The  boundary  conditions  in  the 
xt-plane  corresponding  to  those  given  by  equations  ( 17)  must  be  satisfied. 
Furthermore,  additional  accuracy  can  be  obtained  if  the  aj’s  are  selected 

so  that  the  approximate  solution  given  by  this  method  has  the  same  value 
of  lower  derivatives  at  y  =  0  and  y  =  5  as  an  exact  solution  to  the 
partial  differential  equations  would  have.  These  values  may  be  obtained 
by  differentiating  equation  (l)  one  or  more  times  with  respect  to  y 
or  t.  For  completeness  the  useful  boundary  conditions  in  the  x,t  coordi 
nates  will  now  be  listed  for  flows  having  zero  ax^al  pressure  gradient 
with  or  without  heat  transfer  or  having  an  axial  pressure  gradient  but 
no  heat  transfer: 

At  t  =  0:  u/Ui  =0  ? 

(Ti/T^)  (T-l/tJ  (a2/ hr2)  (u/uj  =  -Pliy&t2 

=0 

At  t  =  6t:  u/ux  =  I 

(^/^(u/uj  =  0(m  =  1,2,3,  .  .  .) 

It  may  be  pointed  out  that  although  additional  conditions  at  y  =  0 
can  be  derived  they  become  inrpractically  cumbersome  to  apply. 
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The  method  involving  the  simultaneous  solution  of  equations  (20) 
and  (21)  will  he  referred  to  as  the  "two -parameter  method/1  while  that 
of  using  only  equation  (20)  will  he  the  "one -parameter  method."  In  the 
two-parameter  method  one  coefficient  a^  in  addition  to  5^  is  deter¬ 
mined  from  the  two  differential  equations,  the  rest  from  the  boundary 
conditions.  In  the  one -parameter  method  h-j.  is  determined  from  the 

solution  of  equation  (20)  and  all  of  the  aj  coefficients  are  deter¬ 
mined  from  boundary  conditions. 

In  carrying  out  the  integration  of  these  differential  equations 
graphical  or  numerical  methods  are  in  general  required.  The  arbitrary 
constants  are  usually  determined  so  that  the  unknowns  are  either  well- 
behaved  or  have  definite  prescribed  values  at  some  value  of  x. 

.  Before  applying  the  methods  indicated  here  to  the  case  of  an  axial 
pressure  gradient,  it  is  considered  desirable  to  investigate  critically 
the  accuracy  obtained  by  using  several  types  of  profiles  in  equation  (20) 
and  in  both  equations  (20)  and  (21)  for  the  compressible  flow  with  no 
axial  pressure  gradient.  An  exact  solution  for  this  case  has  been  given 
by  Chapman  and  Rube  sin  (reference  21).  The  accuracy  of  the  results  will 
be  measured  by  two  criteria,  namely,  by  the  accuracy  of  the  skin-friction 
and  heat-transfer  values  and  by  the  critical  Reynolds  number  for  the 
stability  of  the  laminar  boundary  layer. 

'  The  application  of  these  integral  methods  to  the  flat  plate  and 
the  results  for  skin  friction  and  heat  transfer  will  now  be  discussed 
in  detail. 


BOUNDARY  LAYER  WITH  ZERO  AXIAL  PRESSURE  GRADIENT 


The  one -parameter  and  two -parameter  methods  described  in  the  pre¬ 
ceding  section  are  here  explicitly  applied  to  the  compressible  flow  over 
a  surface  with  zero  axial  pressure  gradient.  Various  profiles  will  be 
assumed  and  the  results  compared  with  those  obtainable  by  a  mathematically 
exact  solution  for  this  case. 


Solution  of  Equations 

For  the  case  of  zero  axial  pressure  gradient  Uj  =  and  equa¬ 
tions  (20)  and  (21)  become 

(8t2)'  =  (2v„CF3)/(Vl) 


(25) 
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and 


(6t2) '  =  (^v_CF6)  /  (26) 

where  the  Fj  integrals  are  in  this  case  constants,  and  not  functions 
of  x. 

In  the  two -parameter  method  equations  (25)  and  (26)  are  solved  for 
and  the  additional  a^  coefficient  which  has  not  been  determined 

from  boundary  conditions.  With  the  initial  condition  5t  =  0  at  x  =  0 
these  unknowns  are  determined  by  the  equations 

2F6/F4  =  F3/Fl  (27) 

and 

8t /*  =  x/(aF3/Fl)  (C/Ex)  (28) 

where  Rx  =  U^x/v^.  Equation  (27)  leads  to  a  quartic  equation  in  the 
unknown  coefficient  a.^,  three  roots  of  which  must  he  rejected  on  physical 
grounds.  Once  the  constant  value  of  a^^  is  known,  F^  and  Fx  can 
he  evaluated  and  substituted  into  equation  (28)  to  complete  the  solution. 

In  the  one-parameter  method  the  Fj  integrals  are  completely  deter¬ 
mined  once  the  velocity  profile  is  chosen  so  as  to  satisfy  the  selected 
boundary  conditions.  Thus,  only  equation  (28)  is  used  to  determine  the 
single  unknown  5^. 


Skin  Friction  and  Heat  Transfer 

From  the  solutions  outlined  in  the  preceding  section,  the  skin- 
friction  and  heat-transfer  characteristics  are  readily  determined.  The 
coefficient  of  average  skin-friction  drag  Cf  is  defined  by  the  equation 


(29) 
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By  use  of  equation  (28)  this  can  be  written  as 

Of  -  2(51*3 (CK)] 1/2  (30) 

Similarly  a  local  heat-transfer  coefficient  h  can  be  defined  as 

h  =  -(k  dT/8y)0|(To  -  Te)  (31) 

where  Te  is  the  equilibrium  wall  temperature  for  no  heat  transfer  and, 
for  Npr  =  1,  is  given  by  the  equation 

(32) 

By  the  use  of  equation  (28)  the  equation  for  h  is  found  to  be 


%u  = 

=  (fxf  3^0/2)  1A  (33) 


Hence,  by  comparison  of  equations  (30)  and  (33)> 

%u-o  A/1*  (3M 

l  ' 

Velocity  Profiles  and  Results  for  Skin  Friction  and  Heat  Transfer 

Several  velocity  profiles  were  assumed  and  the  friction  and  heat 
transfer  given  hy  these  profiles  were  calculated.  In  the  two-parameter 
method  fourth-  and  fifth-degree  polynomials  were  used,  while  in  the  one- 
parameter  method  fourth-,  fifth-,  and  sixth-degree  polynomials  were 
chosen. 

The  results  obtained  here  are  compared  mainly  with  those  of  the 
exact  analysis  of  Chapman  and  Rube sin  (reference  21)  which  treats  a 
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general  Prandtl  number  and  variable  surface  temperature.  The  exact 
velocity  profiles  of  reference  21  are  derived  from  the  well-known  Blasius 
differential  equation,  solutions  of  which  are  tabulated,  for  example, 
in  reference  22*  With  the  assumptions  Npr  =  1  and  constant  surface 

temperature  for  which  equations  (9)  and  (10)  are  valid,  exact  temperature 
profiles  are  thus  obtainable  from  the  exact  velocity  profiles  of  refer¬ 
ence  21.  (These  profiles,  when  expressed  in  terms  of  the  Blasius  vari¬ 
able  T)  (cf.  equation  lk ,  reference  21),  are  independent  of  the  surface 
temperature  and  the  Prandtl  number.) 

The  variation  of  skin  friction  with  Mach  number  for  an  insulated 
surface  is  represented  in  figure  1.  Included  for  comparison  are  results 
already  calculated  in  reference  21  for  Npr  =  0.72,  as  well  as  exact 

results  for  Npr  =  1. 

Table  I  gives  values  of  the  skin-friction  coefficient  and  Nusselt 
number  obtained  by  the  one-  and  two  •-parameter  methods  with  the  various 
velocity  profiles.  Values  derived  from  the  exact  analysis  of  refer¬ 
ence  21  for  Npr  =  1  and  0.72  are  listed  for  comparison.  The  free- 

stream  temperature  and  the  wall  temperature,  which  for  an  insulated 
surface  is  a  function  of  the  free-stream  Mach  number,  are  contained 
explicitly  in  the  term  C.  Typical  values  of  Vc"  for  an  insulated 
surface  are  presented  in  table  II  in  connection  with  the  calculation  of 
the  skin-friction  drag  coefficient. 

Table  I  shows  that  the  values  obtained  for  skin  friction  and  heat 
transfer  by  the  two-parameter  method  with  either  of  the  two  profiles 
used  and  by  the  one -parameter  method  with  a  sixth-degree  profile  differ 
from  the  exact  values  by  less  than  0.6  percent.  The  values  obtained  by 
the  one -parameter  method  with  the  fourth-degree  profiles  are  in  error 
by  approximately  3  percent.  Thus,  all  methods  used  here  give  sufficiently 
accurate  results  for  engineering  purposes  for  skin-friction  and  heat- 
transfer  coefficients. 

Since  in  both  the  exact  solutions  (reference  21)  and  the  solutions 
obtained  here  by  the  integral  methods,  the  effect  of  Mach  number  in  the 
case  of  an  insulated  plate  with  Npr  =  1  is  given  by  the  same  factor, 

namely  \fc,  it  follows  that  (for  a  linear  temperature-viscosity  relation) 
these  integral  methods  lead  to  the  exact  effect  of  Mach  number  on  Cp 

and  %u.  It  may  be  noted  that  the  Mach  number  effect  depends  solely  on 

the  assumed  temperature -viscosity  relation.  The  effect  of  the  latter 
can  be  seen  from  figure  1,  where  the  results  of  the  various  integral 
methods  and  of  reference  21,  all  based  on  a  linear  temperature -viscosity 
relation,  are  shown  along  with  Crocco's  results  which  are  based  on  the 
Sutherland  formula. 
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From  figure  1  it  can  be  seen,  incidentally,  that  the  value  of  the 
Prandtl  number  has  a  small  effect  on  the  variation  of  with  Mach 

number  in  the  case  of  zero  heat  transfer.  For  example,  in  the  range 
M  =  0  to  M  =  with  =  648°  R,  Cf  \fK^  varies  from  1.328  to  O.983 
for  Npr  =  0.72,  and  from  1.328  to  0.954  for  Npr  =  1.  It  can  be  seen 

from  table  I,  moreover,  that  the  Prandtl  number  has  a  somewhat  larger 
effect  on  the  value  of  Nusselt  number  at  all  Mach  numbers,  the  change 
from  Npr  =  1  to  Npr  =0.72  being  about  10  percent. 


Comparison  of  Velocity  Profiles  and  Derivatives  in  Physical  Plane 

For  estimating  the  stability  of  the  laminar  boundary  layer,  the 
accuracy  of  the  first  and  second  derivatives  of  the  velocity  profile, 
as  well  as  that  of  the  profile  itself,  are  of  importance.  In  this  sec¬ 
tion  the  profiles  and  derivatives  obtained  by  the  integral  methods  under 
discussion  are  compared  with  the  exact  Blasius  solution  of  references  21 
and  22.  The  comparison  can  be  conveniently  made  in  the  physical  xy-plane 
by  introduction  of  the  dimensionless  variable 


e  =  (y/^^Rx/c)1/2 


(35) 


The  derivatives  with  respect  to  €  are  shown  in  figures  3  and  4. 

It  is  seen  from  figures  2,  3>  and  4  that,  although  the  velocity 
profiles  obtained  by  the  integral  methods  agree  well  with  the  exact 
profile,  the  derivatives,  as  might  be  expected,  do  not  correspond  so 
closely.  Previous  comparison  of  skin  friction  and  heat  transfer,  which 
depend  upon  the  first  derivatives  at  the  surface,  indicates  that  each 
of  the  integral  methods  gives  good  results  for  these  characteristics. 
Hence  each  of  the  integral  methods  may  be  considered  to  predict  the 
value  of  the  first  derivative  at  the  surface  with  sufficient  accuracy. 
Moreover,  except  for  the  fourth-degree  profiles,  the  first  and  second 
derivatives  given  by  all  the  approximate  solutions  are  in  good  quali¬ 
tative  agreement  with  the  exact  results  throughout  the  boundary-layer 
thickness.  Examination  of  figures  3  and  4  for  first  and  second  deriva¬ 
tives  over  the  boundary -layer  thickness  indicates,  however,  that  the 
one-parameter  method  with  a  sixth-degree  velocity  profile  gives  on  the 
whole  the  most  satisfactory  results.  It  is  to  be  anticipated  that  this 
method  will  consequently  also  give  satisfactory  results  for  stability 
calculations. 
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DETERMINATION  OF  CRITICAL  REYNOLDS  NUMBER  FOR  STABILITY 


The  experimental  work  of  Schubauer  and  Skramstad  (reference  23)  has 
clearly  established  that  the  transition  from  a  laminar  to  a  turbulent 
boundary  layer  is  due  to  an  instability  of  the  laminar  layer,  if  the 
turbulence  in  the  free  stream  is  low,  and  if  the  surface  of  the  body 
is  aero dynamic ally  smooth  and  has  a  large  radius  of  curvature.  Thus 
the  laminar-boundary-layer  stability  theory,  which  has  been  developed 
for  incompressible  flow  by  several  investigators  over  the  past  25  years 
and  which  has  been  recently  extended  by  Lin  and  Lees  (references  4  and  5) 
to  the  practically  interesting  case  of  compressible  flow,  may  be  used 
to  predict  the  local  Reynolds  number  at  which  instability  of  the  laminar 
layer  will  first  occur,  or  above  which  self -propagated  disturbances  are 
not  damped  out.  Transition  to  turbulent  flow  will  take  place  downstream 
of  the  point  corresponding  to  this  local  Reynolds  number,  which  is  termed 
the  critical  Reynolds  number;  the  exact  distance  downstream  cannot  be 
predicted  by  the  present  stability  theory,  which  is  based  on  small- 
perturbation  methods,  but  seems  to  be  dependent  on  the  value  of  the 
critical  Reynolds  number,  on  the  magnitude  of  the  small  but  possibly 
finite  turbulence  in  the  free  stream,  on  the  surface  roughness,  and  on 
the  potential -flow  pressure  gradient. 

The  exact  calculation  of  the  critical  Reynolds  number  for  a  given 
velocity  profile  is  tedious.  However,  the  approximate  stability  rules 
of  Lin  and  Lees  (references  4  and  5)  permit  a  rapid  determination  of 
this  critical  value  with  a  minimum  of  labor.  These  rules  have  been 
shown  by  the  work  of  Lin  and  Lees  and  that  of  Hahneman,  Freeman,  and 
Finston  (reference  24)  to  give  reliable  results  and  have,  therefore, 
been  used  extensively.  It  might  be  further  mentioned  that  because  of 
approximations  inherent  in  the  laminar -boundary-layer  stability  theory, 
and  because  of  the  undetermined  relation  between  boundary-layer  neutral 
stability  and  transition,  there  seems  to  be  no  practical  justification 
for  obtaining  greater  accuracy  than  that  given  by  the  approximate  sta¬ 
bility  rules. 

The  results  of  stability  calculations  are  sensitive  to  the  accuracy 
of  the  profiles.  Thus  the  critical  Reynolds  numbers  based  on  approximate 
solutions  to  the  boundary-layer  equations  have  not  in  general  been  in 
good  agreement  with  those  based  on  exact  solutions.  The  reason  for  this 
discrepancy  is  clear  when  it  is  considered  that  the  stability  calculation 
depends  strongly  on  first  and  second  derivatives  of  the  velocity  profile. 
Approximate  solutions  may  be  expected  to  give  fairly  accurate  velocity 
profiles,  but  unless  special  care  is  exercised  the  first  and  second 
derivatives  of  the  profile  throughout  the  entire  boundary-layer  thick¬ 
ness  might  be  quite  in  error.  This  has  been  shown  in  the  previous  sec¬ 
tion  where  it  has  been  pointed  out  that  the  frequently  used  Karman- 
Pohlhausen  method  with  a  fourth-degree  velocity  profile  gives  reasonably 
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accurate  values  for  the  skin  friction,  the  heat  transfer,  and  the  velocity 
profile  hut  yields  quite  inaccurate  values  for  the  derivatives.  It  may  m 

he  -expected,  therefore,  that  a  crucial  test  of  the  accuracy  of  approxi¬ 
mate  solutions  to  the  laminar -boundary -layer  equations  would  he  a  com¬ 
parison  of  the  critical  Reynolds  numbers  predicted  by  these  solutions 
with  those  predicted  by  exact  solutions. 

In  this  section  the  critical  Reynolds  numbers  of  the  profiles 
obtained  in  the  previous  section  for  the  zero-heat-transfer  case  are 
obtained  and  compared  with  those  of  the  exact  solution  of  Chapman  and 
Rube  sin. 

With  the  prime  denoting  partial  differentiation  with  respect  to  e, 

Lees*  criterion  modified  for  the  viscosity-temperature  relation  used  here 
is  given  as  follows: 


Rcr  ~  2(Rx,crc)1/^2 

-  25(u/°l)o,C(T/Tl)2(U/Ul)'1‘{1  -  Ml2£  -  ("/Ul)^j  2}  -1/2  W 


where  the  subscript  c  denotes  values  at  the  point  6  =  ec  at  which 
the  following  relations  are  satisfied: 

J  -  -”(Tl/To)(u/Ul)o'(U/Ul)[(',/Ul)]'3[(T/Tl)("/Ul)''  -  (”/%)  '(T/*l)] 

(37) 

and 

0.580/J  =  1  -  2|{u/U1)o,€(u/U1)  -  lj  (38) 

It  will  be  noted  from  equation  (37)  that  the  first  and  second  deriva¬ 
tives  of  the  velocity  and  temperature  profiles  in  the  critical  region 
around  ec  are  influential  in  the  determination  of  Rcr.  The  procedure 

for  obtaining  Rcr  is  to  calculate  J  from  equation  (37)  for  various 
values  of  e,  and  to  find  the  values  of  €  =  ec  for  which  equation  (38) 
is  satisfied.  After  the  values  of  (u /^i)c  and  (T/Tl)c  are  determined 
at  this  point,  Rcr  and  Rjj  cr  follow  from  equation  (36). 
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For  comparison  with  previous  results  it  is  convenient  to  calculate 
the  critical  Reynolds  number  R0  cr  based  on  momentum  thickness.  Since 

0  is  defined  by 


(p/pi)W“i)[1  -  W“0] 


(39) 


it  is  found,  for  the  integral  methods,  that 

Re,cr  =  (FlF3/2)l/2Rcr  (^) 


and  for  the  profiles  of  reference  (21) 


R0,cr  =  0-332Rcr 


m 


Calculations  of  the  minimum  critical  Reynolds  numbers  have  been 
made  for  the  flow  over  a  flat  insulated  surface  at  free -stream  Mach 
numbers  of  0,  1,  and  2.  The  results  are  tabulated  in  table  III  and 
plotted  in  figure  Primarily,  comparison  of  the  integral  methods 

presented  here  should  be  made  with  the  results  derived  from  the  exact 
profiles  of  reference  21  since  identical  temperature -viscosity  relations 
are  used  in  both.  The  results  of  Lees,  given  in  figure  5  of  reference  5> 
are  also  presented  for  comparison.  Note  should  be  taken  that  Lees* 
approximate  values  are  based  on  still  a  further  approximation  to  equa¬ 
tions  (36),  (37)>  and  (38)  presented  here. 

Figure  5  indicates  that  the  one-parameter  sixth-degree  polynomial 
method  gives  the  best  general  agreement  with  results  of  the  exact  solu¬ 
tions  over  the  range  of  Mach  numbers  studied.  The  two -parameter,  fifth- 
degree  method  is  next  best,  while  in  order  of  decreasing  accuracy,  the 
one -parameter  fifth-degree,  two -parameter  fourth-degree,  and  one -parameter 
fourth-degree  polynomials  follow.  The  two -parameter  methods  seem  to 
give  better  general  agreement  than  the  one -parameter  methods  for  equal 
degree  of  the  assumed  polynomial  for  the  velocity  profile.  (Because  of 
the  inaccuracies  introduced  by  the  approximation  in  this  stability 
criterion,  as  well  as  by  the  required  graphical  method  involved  in  its 
application,  the  derived  results  are  not  reliable  to  more  than  two  signif¬ 
icant  figures.) 
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The  stability  analysis  given  here  has  thus  indicated  that  integral 
methods  of  solution  of  the  boundary-layer  equations  may  be  used  to 
predict  the  critical  Reynolds  number  provided  that  either  additional 
integral  equations  or  higher-degree  polynomials  and  additional  boundary 
conditions  with  a  single  integral  equation  are  used.  Thus  the  choice 
of  method  to  be  used  in  the  more  general  cases  of  pressure  gradient  with 
and  without  heat  transfer  is  either  the  two-parameter  method  or  the  one- 
parameter  method  with  higher-degree  polynomials,  particularly  of  the 
sixth  degree. 

While  these  two  methods  lead  to  results  which  are  of  nearly  the 
same  accuracy,  use  of  the  one -parameter  method  has  the  advantage  of 
simplicity  of  calculation,  which  becomes  significant  in  cases  involving 
pressure  gradient.  It  is  thus  concluded  that  the  one -parameter  method 
with  sixth-  or  higher-degree  polynomials  is  the  most  satisfactory  and 
promising  method  for  extension  to  general  compressible  boundary-layer 
analyses.  Sixth-degree  polynomials  will  probably  prove  to  be  satis¬ 
factory  in  many  cases,  but  use  of  a  seventh-degree  profile  satisfying 
an  additional  condition  at  the  separation  point  (cf.  reference  9)  may 
increase  the  accuracy  of  determining  the  separation  point  in  adverse- 
pressure  -gradient  cases.  Moreover,  a  sixth-degree  profile  may  not  be 
satisfactory  for  flow  near  a  stagnation  point  (cf.  reference  7)« 


CONCLUSIONS 


The  laminar -boundary-layer  equations  for  compressible  flow  can  be 
converted  into  one,  two,  or  more  integral-differential  equations.  Approxi 
mate  solutions  qan  be  obtained  by  assuming  special  forms  for  the  velocity 
profiles  satisfying  various  boundary  conditions  and  containing,  in 
practice,  either  one  or  two  parameters  to  be  determined  by  these  equations 

From  a  comparison  of  the  integral  methods  discussed  here  for  the 
compressible  flow  over  a  surface  with  no  axial  pressure  gradient,  the 
following  conclusions  can  be  drawn: 

1.  All  of  the  methods  predict  the  values  of  the  skin-friction  and 
heat-transfer  coefficients  at  the  wall  as  well  as  the  over-all  velocity 
profiles  with  satisfactory  accuracy  for  engineering  purposes.  The  two- 
parameter  method  and  the  one -parameter  method  with  a  sixth-degree  pro¬ 
file  are  particularly  accurate. 

2.  The  one -parameter  method  with  a  sixth-degree  velocity  profile 
gives  on  the  whole  the  most  accurate  results  for  velocity  profiles  and 
their  derivatives  throughout  the  boundary-layer  region. 
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3.  The  critical  Reynolds  number  for  laminar -boundary-layer  sta¬ 
bility  is  predicted  with  qualitative  accuracy  by  all  methods.  Moreover, 
reasonably  good  quantitative  accuracy  is  obtained  with  the  two -parameter 
method  and  with  the  one -parameter  method  with  a  sixth-degree  profile. 

Because  of  computational  simplicity  and  equality  of  accuracy, 
the  one -parameter  method  with  sixth-  or  higher-degree  polynomials  appears 
to  be  the  most  satisfactory  and  promising  method  of  investigating  the 
more  general  case  of  laminar  compressible  boundary  layer  with  axial 
pressure  gradient  and  heat  transfer. 


Polytechnic  Institute  of  Brooklyn 

New  York,  N.  Y.,  November  6,  1950 
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TABLE  I.-  SKIN -FRICTION  DRAG  COEFFICIENT 
AND  NUSSELT  NUMBER 


Method 

Degree  of 
polynomial 

Boundary  conditions 
satisfied 

(cf.  equations  (24)) 

1  /M1/2 

Hi) 

2NNu(RxC)1/2 

Two  parameter 

4 

0.663 

0.663 

5 

.660 

.660 

u 

l,2,4,5(m  =  1,2) 

0.685 

0.685 

One  parameter 

5 

l,2,3,4,5(m  =  1,2) 

.644 

.644 

6 

1,2,3, 4, 5(m  =  1,2,3) 

.661 

.66l 

Reference  21 

TQ  =  Constant 
NPr  =  1 

0.664 

0.664 

Reference  21 

T0  =  Constant 
NPr  =  0.72 

____ 

0.664 

0.592 

- - - - - - - 
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TABLE  II.-  VALUES  OF  \fc  FOR  AN 


INSULATED  SURFACE 

FOR  Npr  =  1 

Moo 

T 

00 

72°  F  absolute 

• 

648°  F  absolute 

0 

1.0 

1.0 

1 

1.021 

.976 

2 

1.057 

.916 

3 

1.074 

-3- 

00 

• 

4 

1.067 

.776 

5 

1.043 

co 

1 — 1 

r- 

• 

example,  2  parameter-4  means  two-parameter  method  with  a  fourth-degree  profile. 
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Figure  1.-  Variation  of  skin-friction  drag  coefficient  with  free-stream 
Mach  number  for  insulated  flat  surface.  I,  one -parameter  method) 

II,  two-parameter  method)  4,  5,  6,  degree  of  polynominals)  C-R, 
Chapman-Rube sin  (reference  21) )  C-S,  Crocco's  calculation  based  on 
Sutherland's  temperature -viscosity  relation  (cf.  reference  21). 
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(a)  M*  =  0. 

Figure  2.-  Comparison  of  velocity  profiles  for  flow  over  insulated  flat 

surface . 
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U, 


0>)  Mqo  =  2. 


Figure  2.-  Concluded. 
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Figure  3.-  Comparison  of  first  derivatives  of  velocity  profiles  for  flow 

over  insulated  flat  surface. 
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Figure  3--  Continued. 
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(c)  Two-parameter  method:  M  =  2. 

^  00 


Figure  3.-  Continued. 


e  =  -X-  /i 

ZXs!  C 


(d)  One-parameter  method;  =  2. 


Figure  3.-  Concluded 
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Figure  -  Comparison  of  second  derivatives  of  velocity  profiles  for  flow 

over  insulated  flat  surface. 
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(c)  Two-parameter  method; 

Figure  4.-  Continued. 
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(d)  One-parameter  method; 

Figure  4.-  Concluded. 
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Figure  5.-  Variation  of  minimum  critical  Reynolds  number  with  Mach  number 
for  laminar  boundary  layer  flow  over  insulated  flat  surface  with 

Npr  =  1. 
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stability.  From  this  comparison,  it  is  found  that  stability.  From  this  comparison,  it  is  found  that 

Karman-Pohlhausen  method  with  the  sixth-degree  Karman-Pohlhausen  method  with  the  sixth-degree 

polynomial  as  the  velocity  profile  is  the  most  suitable  polynomial  as  the  velocity  profile  is  the  most  suitable 

for  many  practical  purposes.  for  many  practical  purposes. 


